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ABSTRACT 

We present hydro dynamic simulations of high- velocity clouds (HVCs) traveling through the hot, 
tenuous medium in the Galactic halo. A suite of models was created using the FLASH hydrodynamics 
code, sampling various cloud sizes, densities, and velocities. In all cases, the cloud-halo interaction 
ablates material from the clouds. The ablated material falls behind the clouds, where it mixes with 
the ambient medium to produce intermediate-temperature gas, some of which radiatively cools to less 
than 10,000 K. Using a non-equilibrium ionization (NEI) algorithm, we track the ionization levels of 
carbon, nitrogen, and oxygen in the gas throughout the simulation period. We present observation- 
related predictions, including the expected H I and high ion (C IV, N V, and O VI) column densities on 
sight lines through the clouds as functions of evolutionary time and off-center distance. The predicted 
column densities overlap those observed for Complex C. The observations are best matched by clouds 
that have interacted with the Galactic environment for tens to hundreds of megayears. Given the 
large distances across which the clouds would travel during such time, our results are consistent with 
Complex C having an extragalactic origin. The destruction of HVCs is also of interest; the smallest 
cloud (initial mass ^120 M ) lost most of its mass during the simulation period (60 Myr), while 
the largest cloud (initial mass %4x 10 5 M ) remained largely intact, although deformed, during its 
simulation period (240 Myr). 

Subject headings: Galaxy: halo — hydrodynamics — ISM: clouds — methods: numerical — turbu- 
lence — ultraviolet: ISM 



1. INTRODUCTION 



The halo of the Milky Way contains clouds of neu- 



tral hydrogen (H I) with |^lsr| > 90 km s~", 
known as high- velocity c l ouds (HVCs; iMuller et al.l 1963: 
Wak ker fc van Woerdenl [1997) . Although some HVCs 
are relatively close to the Galactic disk (e.g., \z\ < -4 kp c 
for Complex M; iDanlv et al.l 119931 : iKeenan et al1ll995D . 
other HVCs are known to be in the upper Galactic halo, 
where the Galaxy interacts wit h its surroundings (e.g., 
d = 10 ±2.5 kpc for Complex C: lThom et all[2008h . The 
distant HVCs are thought to be gas stripped off satellite 
galax ies (e.g., iGar diner fc Noguchil 119961 : iPutman et al. 
2004) . extragalactic g as falling into the Galaxy ([Port 
[19661 : IMit^t^[T999h. or gas left over fro m the forma- 
tion of the Galaxy (jMaller fc Bullock|[200l . 

Observations of interstellar absorption of light from 
AGN show that the halo also contains highly ionized 
gas (including high ions such as C IV, N V, and O VI) 
movin g with velocities comparable to those of the H I 
HVCs (E embach eTHI [20031 : IFox et al.ll200i [20051 . [20061 : 
iCollins et al.l 120071 ). These high ions trace gas with 
T ~ (1-3) x 10 5 K, and enable us to probe the connection 
between the hottest Galactic gas (T > 10 6 K) and cooler 
phases of the interstellar medium (ISM). These ions also 
help us investigate the interaction between HVCs and the 
ambient halo gas, and, by extension, the properties of the 
clouds and the halo gas. Some highly ionized HVCs are 
seen in the same directions as H I HVCs, but there are 
also highly io nized HVCs without corresponding high- 
velocity H I (jSembach et all 12003). The sky covering 
fraction of O VI HVCs f~60%: ISembacli_elalJ[20(^ is 
larger than that of H I HVCs (^37%; [Lockman et al.l 



2002), suggesting the possibility that in the halo, warm- 
hot (a few times 10 5 K) high-velocity gas traced by high 
ions could be more common than cool high- velocity gas 
traced by H I. 

High ions in HVCs could hypothetically result from 
several physical processes, including mixing of cool neu- 
tral HVC gas with hot (T > 10 6 K) highly ionized am- 
bient gas, thermal conduction between cool HVC gas 
and hot ambient gas, radiative cooling of hot high- 
velocity gas, and shock heating. For example, when 
a cool cloud is embedded in hot gas, thermal conduc- 
tion can form a layer of intermediate temperature gas 
(T ~ a fe w times 10 5 K) in whi ch high ions are abun- 
dant (e.g.. iBorkowski et al.lll990l for stationary clouds). 
Also, simple radiative cooling of hot gas will result in it 
passing through the temperature range optimal for high 
ions. Collisional ionization equilibrium (CIE) and non- 
equilibrium ionization (NEI) calculations of this cool- 
ing, for stationary g as at least, have been c arried out 
by many authors (see lGnat fc Sternbergll2007l and refer- 
ences therein). High ions could also be produced in the 
hot gas behind the shocks that form as HVCs travel su- 
personically through the ambient medium. The observed 
ratios of high ions' column densities or line intensities can 
be compared with the results of these models to deter- 
mine how high ions are produced; e.g., see Figure 7 in 
KS10. Although this particular comparison was made for 
low- velocity ions, one could test models for high- velocity 
ions in a similar way. 

Some of these hypothetical processes require 
that hot gas is plentiful in the halo. Observa- 
tions of the diffuse s o ft X-r a y background (S XRB; 
iBurrows fc Mendenhalll 119911 : ISnowden etHI 119911 . 
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1991; Emtz fc Snowdenl 1200(1 [Smith et al.l 120071: 
Galeazzi et al. 2007; Henlev & Shelton 2008; Lei et al. 
20091 : lYoshino et all 120091 : iHenlev et all l2010f ) and of 
X-ray absorption lines, su c h as O VII and O VIII 
(e.g., lYaofcWand [2QQpL Ig O DTt iFang et aH 120061 : 
iBregman fc Llovd-Daviesl 120071 : lYao et al.l I2009D show 
that hot gas with T ~ 1-3 x 10 6 K is indeed present in 
the halo. Although the distance and origin of this hot 
gas are uncertain at present, it is likely that some of this 
hot gas exists in the same region as the H I HVCs. If so, 
then it is plausible that high ions are produced in the 
turbulent mixing layers (TMLs) between the cool HVCs 
and the hot ambient ISM. This turbulent mixing of cool 
and hot gas arises from the Kelvin-Helmholtz or shear 
instability induced by the velocity difference across the 
interface between the t wo types of gas. The TML m odel 
was first suggeste d by [B egelman & Fabian (1990) and 
later developed by Sl avin et al . (1993). Column densities 
of various high ions (including C IV, N V, and O VI) 
have been predicted for the TML model us ing a variety 
of te chniques: analytical calculations (|Slavin et al.l 
1993) , 3D magnetoh ydrodynamic (MHD) simulations 
(jEsquivel et al.l 12006). and 2D hydro dynamic simula- 
tions incorporating NEI calculations (jKwak fc Sheltor] 
I2010L hereafter KS10). Note that these studies all 
used a plane-parallel geometry, rather than a HVC-like 
geometry. TML column density predictio ns have been 
compared with obser v ations of HVCs (jSembach et al.l 
12003 IFox et all l200l iCollins et all 120071 : KS10). and 
should be of value for studying other astrophysical 
situations in which cool or warm material slides past 
hot material. 

In this paper, we present the results of simulations 
that trace the interaction between an H I HVC and the 
hot ambient medium in a more realistic geometry than 
previous TML studies, that is, a spherical cloud falling 
through the hot ISM. As in our previous study (KS10), 
we trace the ionization states of carbon, nitrogen, and 
oxygen with non-equilibrium calculations which allow us 
to estimate the amounts of interesting ions (C IV, N V 
and O VI) more accurately than CIE calculations. The 
new simulations enable us to investigate how the cold 
gas ablates from a spherically shaped cloud due to shear 
instabilities, how it mixes with the hot ambient gas, and 
where the high ions that are produced in the process of 
ablation and mixing reside. In addition, we can examine 
the velocities of the high ions: there are high- velocity 
high ions that move with similar velocities as the H I 
HVCs and low- velocity high-ions that move with similar 
velocities as the ISM. Using the results of our simula- 
tions, we will answer the following questions: (1) How 
many high ions are produced by H I HVCs traveling 
through hot ambient ISM? (2) How many high-velocity 
high ions result? (3) How do the ratios of high ions and 
H I compare with observations? (4) Are HVCs likely 
to reach the Galactic disk intact? In this paper we will 
concentrate on the high- velocity ions, deferring the dis- 
cussion of low- velocity ions to Paper II (D. B. Henley 
et al. 2011, in preparation). 

This paper is organized as follows. In the next section, 
we briefly summarize our numerical methods and present 
the parameters for our suite of seven simulational mod- 
els. In Section[3]we present the results of our simulations. 



Specifically, Section 13.11 describes in detail the hydrody- 
namical evolution of one of our model clouds, Section 1X21 
describes the evolution of the amount of neutral mate- 
rial, Section [3~3l describes the effect of the di ffere nt model 
parameters on the cloud evolution, Section 13.41 discusses 
whether or not HVCs can reach the Galactic disk, and 
Section 13.51 describes the production of high ions in our 
simulations. In Section [4] we compare the column den- 
sities and column density ratios predicted by our simu- 
lations with observed Complex C values. In Section [5] 
we discuss the effects that the various assumptions in 
our simulations have on our results. We summarize our 
results in Section [6] 

2. NUMERICAL METHODS AND MODEL PARAMETERS 

We carried out our simulations using the same code a s 
KS10, namely FLASH version 2.5 (jFrvxell et al.ll2000l ). 
We used the FLASH NEI module to track the ionization 
evolution of carbon, nitrogen, and oxygenQ The sim- 
ulations include radiative cooling, although the cooling 
curve was calculated assuming CIE. See KS10 for more 
details of the code, including some discussion of CIE ver- 
sus NEI cool ing ra t es. A s in KS10, we used the abun- 
dances from Allen (1973), which are the default abun- 
dances in FLASH. Our results can be rescaled to give 
results for different abundance tables, assuming that the 
radiative cooling curve does not change significantly with 
the assumed abundances (an assumption which may not 
always be valid). 

In KS10 the simulations were run in 2D Cartesian co- 
ordinates to study TMLs in a plane-parallel geometry. 
Here, we use 2D cylindrical coordinates to study initially 
spherical clouds as they fall through the hot ISM. As 
in KS10, we use 2D simulations in order to minimize 
the unaffordably large computing resources (particularly 
memory) needed to track the ionization of the three el- 
ements in our NEI simulations. In addition, instead of 
the cloud moving relative to a stationary ISM, in our 
simulations the ambient medium moves relative to an 
initially stationary cloud. In this way, we can trace the 
cloud evolution for a long time without requiring a large 
computational domain, reducing the amount of memory 
and CPU time needed to run the simulations. However, 
although the cloud is initially stationary in the compu- 
tational domain, all velocities stated in this paper are 
in the observer's frame, in which the ISM is stationary. 
In the computational grid's reference frame, the ISM is 
initially moving in the +z direction with velocity |i> ZjC i|, 
where — \v ZiC \\ is the initial velocity of the cloud in the ob- 
server's frame. Therefore, the z- velocity of gas in a given 
grid cell in the observer's frame is v z = v z , s i m — |^z,ci|, 
where v z , s i m is the z- velocity of the gas in that cell ob- 
tained directly from the simulation. Note that, although 
we present results for HVCs with negative velocities (i.e., 
moving toward us), our results are equally applicable to 

1 Si IV is another high ion that has been observed in the halo 
of the Milky Way via its UV lines. However, this ion is more 
susceptible to photoionization than C IV, N V, or O VI, due to 
its lower ionization potential. Modeling photoionization is beyond 
the scope of this paper, so we do not include silicon in our NEI 
calculations. Note that the ionization and recombination rates in 
the FLASH NEI module include only the effects of collisional ion- 
ization, auto-ionization, radiative recombination, and dielectronic 
recombination. 
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Table 1 

Model Parameters 



Model 
(1) 


Radius 
(pc) 
(2) 


V z ,cl 

(km s- 1 ) 
(3) 


™H,cl 

(cm- 3 ) 
(4) 


M init ,T 

(M ) 
(5) 


-^init,v 

(M ) 
(6) 


™H,ISM 

(cm" 3 ) 
(7) 


A 


20 


-100 


0.1 


120 


130 


1.0 x 10- 4 


B 


150 


-100 


0.1 


4.9 x 10 4 


5.1 x 10 4 


1.0 x 10- 4 


C 


150 


-150 


0.1 


4.9 x 10 4 


5.1 x 10 4 


1.0 x 10- 4 


D 


150 


-300 


0.1 


4.9 x 10 4 


5.1 x 10 4 


1.0 x 10- 4 


E 


150 


-150 


0.1 


4.9 x 10 4 


4.9 x 10 4 


1.0 x 10- 4 


F 


300 


-100 


0.1 


4.0 x 10 5 


4.2 x 10 5 


1.0 x 10- 4 


G 


150 


-100 


0.01 


4.9 x 10 3 


5.1 x 10 3 


1.0 x 10- 5 



Note. — Col. (1): Model identifiers. Col. (2): Approximate radius of the 
model cloud (except for Model E). See Figure [J for the detailed density profile 
of each model cloud. Model E is a uniform-density cloud with an exact radius. 
Col. (3): Initial velocity of the cloud along the z-direction measured in the 
observer's frame. Col. (4): Initial hydrogen number density of the cloud at 
its center. Col. (5): Initial mass of cloud having a temperature T < 10 4 K. 
Col. (6): Initial mass of cloud moving with v z c \. Note that all cloud material 
with a hydrogen number density greater than 5riH,iSM: where tih,ism is the 
hydrogen number density of the ISM (Col. 7), moves initially at speed v z ^ c \ 
relative to the ISM. All models have the same cloud temperatures (T c \ = 10 3 K 
at cloud centers) and ISM temperatures (Tjsm = 10 6 K). 

HVCs with positive velocities if the signs on the velocities 
are changed. 

The parameters of our various model clouds are pre- 
sented in Table Q] and Figure [TJ We ran seven different 
models in which we varied the cloud's initial radius (Ta- 
ble [TJ Col. 2), radial profile (Figure [TJ, z- velocity, v z ^ c \ 
(measured in the observer's frame; Table [TJ Col. 3), and 
number density expressed in terms of the hydrogen num- 
ber density, where nn/^He = 10 (Table [TJ Col. 4). The 
ambient number density is 10 3 times smaller than the 
initial cloud number density (Table [TJ Col. 7). All cloud 
material with a density > 5 times the ambient density 
initially moves at v ZjC \ relative to the ambient medium, 
while the less dense outskirts of the cloud are initially 
set to the ISM's velocity. As in KS10, we required that 
the cloud and the ISM were initially in pressure bal- 
ance. We achieved this by varying the temperature from 
T cY = 10 3 K at the cloud center to Tt S m = 10 6 K in the 
ambient medium, such that the pressure was constant. 

Because the clouds generally do not have sharp edges, 
the cloud mass is not well defined. For this reason, we 
have calculated two masses for each cloud: the mass of 
material with T < 10 4 K (Table HJ Col. 5), and the mass 
of material initially moving at v z ^ c \ relative to the am- 
bient medium (Table [TJ Col. 6). These two masses are 
generally similar to each other. The cloud in Model E 
has a sharp edge, and so the two masses are identical in 
this model. 

By running Models B, C, and D, we can see the effects 
of varying the cloud's initial velocity (Section l3.3.ip . be- 
cause these models otherwise have the same initial pa- 
rameters. Apart from the radial density profile, Models C 
and E are identical, and so these two models sh ow th e 
effect of the cloud's initial density profile (Section [3. 3. 2|) . 
The cloud and ISM densities in Model G are an order of 
magnitude smaller than those in the other models, and so 
comparing Models B and G reveals the effect of varying 
the cloud's density (Section 13.3.2)) . Finally, Models A, 
B, a nd F r eveal the effect of varying the cloud size (Sec- 
tion E33J. 




100 200 300 400 500 600 
radius (pc) 

Figure 1. Initial hydrogen number density profiles for our various 
models (see also Table [TJ. 

Although the number densities of HVCs and of the 
upper halo are not well constrained, our chosen model 
parameters are consistent with existing observational 
and theoretical constraints. The H I column densi- 
ties on sightlines passing through the centers of our 
model clouds range from 9.3 x 10 18 cm -2 (Model G) 
to 1.9 x 10 20 cm -2 (Model F), and are consistent 
with measured column d ensities of H I HVCs (e.g., 
iHulsbosch fe Wakkerl I1988D . The density of the hot 
gas in the upper halo is more uncertain, with observa- 
tional estimates including nn < 6.3 x 1Q~ 4 cm" 3 from 
pulsar dispersion measures (GaenslereTaD |2008) and 
Tin — 7.5 x 10 -4 cm -3 from O VII column density mea- 
surements (assuming a uniform spherical halo of radius 
20 kpc; Bregman fc Llovd-Dav ies 2007)0 The halo den- 
sity has also been constra ined by combining observations 
with theoretical models. iPeek et all ((2007) obtained an 
upper halo density of nn ~ 2 x 10 -4 cm -3 by consider- 
ing the drag force on an H VC complex observed in H I. 
iGrcevich fc Putmanl (|2QQ9h estimated that the number 

2 These estimates of the density of the upper halo were given as 
electron number densities. We have converted them to hydrogen 
number densities assuming nn/^He — 10- 
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density in the upper halo is > (2-3) x 10 -4 cm -3 , assum- 
ing that the Milky Way's satellite galaxies lost their gas 
through ram-pressure stripping as they passed through 
the upper halo. The halo density may not decrease dra- 
matically with distance - Wei ner fc Williams! (|1996[ ) es- 
timated that the density near the Magellanic Stream 
(d ~ 50 kpc) is nn ^ 1 x 10 -4 cm -3 , assuming that 
the observed Ha emission arises from an interaction be- 
tween the Magellanic Stream and the ambient gas. Note 
that most of these estimated halo densities are slightly 
larger than the halo density used in most of our simu- 
lations ( nn = 1 x 1Q~ 4 cm" 3 f or Models A through F). 
However, Sembac h et al.l ([20031 ) found that their observa- 
tions of O VI HVCs favored a low-density extended hot 
halo (n H < 10" 4 -10" 5 cm -3 , r > 70 kpc). The lower 
halo density in Model G (nn = 10 -5 cm -3 ) was chosen 
to investigate the lower range of halo densities, as well 
as a lower cloud density. 

3. RESULTS 

3.1. Hydrodynamical Evolution of High- Velocity Clouds 

We chose Model B as our representative model as it 
exhibits many of the physical processes also seen in the 
other models. In Section [3~3], below, we will discuss the 
differences and similarities between the various models. 
Here, we discuss Model B, whose evolution is shown in 
Figure [2j Figure [2] plots, from top to bottom, the hy- 
drogen number density, the temperature, v z (in the ob- 
server's frame), and the ion fraction^) of C IV, N V, and 
O VI. Note that at t = Myr the ion fractions for all 
three high ions are low almost everywhere on the grid. 
This is because the cloud is too cold and the ambient 
medium too hot for these ions. 

Figure [2] shows that the initially spherical cloud de- 
forms during the simulation. This is due to the Bernoulli 
effect. As the ISM flows around the cloud, the speed of 
the ISM relative to the cloud is greater along the sides 
of the cloud (r ~ 150, z ~ pc) than near the top or 
the bottom of the cloud (r ~ 0, z ~ ±150 pc). As a re- 
sult, according to Bernoulli's equation, the pressure will 
be greater immediately in front of and behind the cloud 
than at the edge. Examining the pressures in the output 
from the hydrodynamical code confirms this expectation. 

The cloud also defo rms due to Kelvin- H elmholtz or 
shear instabilities (e.g. JChandrasekharlll961[ ), instigated 
by the velocity difference between the cloud and the ISM. 
Because the speed of the ISM relative to the cloud is 
largest at the edge of the cloud, the ins tabilities grow 
most rapidly here ([Chand rasekharl 1 19611 Section 101), 
and this part of the cloud is pulled outward (see the hy- 
drogen number density and v z plots in Figure [2]) . As the 
ISM flows around the initially spherical cloud, a vortex 
develops behind the cloud. The flow in this vortex be- 
comes more complicated as the shear instabilities start to 
ablate material from the edge of the cloud. This material 
impedes the flow of the ISM into the vortex behind the 
cloud, while new vortices form around the ablating ma- 
terial. This complicated flow causes ISM material some 
way behind the cloud to move in toward the r = axis; 
this material then flows down (i.e., in the — z direction) 

3 For example, the C IV ion fraction is the fraction of all carbon 
atoms that are in the C +3 ionization state. 



and then out again. As this material flows outward along 
the back of the cloud, it helps to stretch the cloud out in 
the horizontal direction. 

The gas ablated from the cool, dense cloud (T c \ = 
10 3 K, 7iH, c i = 0.1 cm -3 ) mixes with the hot, tenuous 
ISM (Tism = 10 6 K, n H ,iSM = 10" 4 cm" 3 ), creating 
mixed gas of intermediate temperature. Eventually, the 
mixed gas reaches a temperature of a few times 10 5 K, 
which is optimal for radiative cooling through line emis- 
sion. As the mixed gas flows back from the edge of the 
cloud, the flow splits: some of the mixed gas is drawn 
into the vortex behind the cloud, while some flows fur- 
ther back from the cloud, continuing to mix with the hot 
ISM. The temperature of the mixed gas flowing away 
from the cloud continues to increase above a few times 
10 5 K, as the continued mixing raises the temperature 
more rapidly than the gas can radiatively cool. In con- 
trast, the mixed gas flowing into the vortex behind the 
cloud cools, due to both mixing with cooler gas and ra- 
diative cooling. In particular, this gas cools more effi- 
ciently after it reaches the region near the r = axis, 
resulting in the accumulation of gas with T ~ 10 4 K 
along this axis in the temperature plots of Figure [2j We 
traced the fractions of original cloud material and ISM 
material contained in the mixed gas and found that a 
significant fraction of the cool gas along the r = axis 
at later times was initially hot ISM, indicating that this 
gas is mixed gas that has undergone radiative cooling. 

The plots of the high ion fractions (the last three rows 
in Figure [2]) show that the fractions of these high ions 
are higher in the mixed gas than in the initially cool 
cloud gas or in the initially hot ISM gas. As the cool 
and hot gas mix, the ions that we are interested in are 
produced both by ionization during the heating of the 
initially cool gas, and by recombination during the cool- 
ing of the initially hot gas. The fraction of Li-like ions 
for each element depends somewhat on the temperature 
of the gas; for example, in the hotter gas the O VI frac- 
tion is higher than the C IV fraction. The NEI ionization 
levels in Figure [2] are different from those expected from 
CIE, because changes in the ionization levels lag behind 
the changes in temperature that are due to mixing or ra- 
diative cooling. This is similar to what we found in our 
previous study of TMLs (KS10). 

It is possible that our use of 2D cylindrical geometry 
resulted in an overestimate of the amount of cool gas 
that accumulates along the r = axis. In such a ge- 
ometry, material that flows toward the r = axis tends 
to stick to the axis, because of the reflecting boundary 
condition there. The 2D cylindrical geometry is also re- 
sponsible for the protuberance at the front of the cloud 
seen at later times, for example, the feature indicated by 
the arrow in the far r ight density plot in our Figure |5] 
(jVieser fc Henslerl [2007bh . However, although the high 
ions contained in the cooled gas that accumulates along 
the symmetry axis give rise to large ion column densities 
along this sightline (see Section [33]) , this column of ma- 
terial does not contribute much to the total number of 
high ions on the grid. 

3.2. Evolution of the Number of H I Atoms 

The default FLASH NEI module assumes that hydro- 
gen is fully ionized, and so we were unable to trace the 
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Figure 2. Cross-sections through the Model B cloud, showing the time evolution of (from top to bottom) the density (expressed in terms 
of the hydrogen number density), temperature, v z (in the observer's frame), C IV ion fraction, N V ion fraction, and O VI ion fraction. 
The hydrogen number density includes both neutral and ionized hydrogen. Each variable is plotted at 15 Myr intervals from t = to 
t = 120 Myr. Note that all variables apart from v z are plotted with logarithmic color scales. The arrow in the far right density plot 
indicates the protuberance mentioned in the text. 



ionization evolution of hydrogen. We therefore assumed 
that the hydrogen in gas with T < 10 4 K is entirely 
neutral (H I), while the hydrogen in hotter gas is fully 
ionized. Our initial clouds "lose" their H I content via 
the physical processes discussed above, namely ablation 
of material from the cloud and its subsequent temper- 
ature increase due to mixing with the hot ISM. (Note, 
however, that H I can also be replenished when hot gas 
cools below 10 4 K, and such gains can count against the 
losses.) Here, we consider the loss of H I from the clouds 
due only to heating, and the loss of H I due to heating 
and/or ablation. In the former case, we consider all H I 
atoms, regardless of their velocity. In the latter case, we 
consider only H I atoms moving with HVC-like velocities: 
v z < -80 km s _1 for Models A, B, F, and G (in which the 
initial velocity of the cloud was v z ^\ = —100 km s _1 ), and 
v z < —100 km s _1 for Models C, D, and E (in which v z ^\ 
was —150 or —300 km s _1 ). These definitions were cho- 
sen to correspond with observational analyses, in which 
one can distinguish between high- and low- velocity H I. 

To investigate the loss of H I from the clouds, for each 
model we calculate the ratio, /?(£), of the number of neu- 
tral hydrogen atoms lost since the beginning of the simu- 
lations to the initial number of neutral hydrogen atoms, 



i.e., 



A/h 



I.init 



■Mn(t) 



■A/h 



I,init 



(i) 



where A/Hi,init and Mm(t) are the initial and current 
numbers of H I atoms, respectively. As noted above, we 
have two different ways of counting A/hi(£), leading to 
two different values of /3(t): /3aii(£) f° r H I at all veloci- 
ties, and /3hvc(^) for HVC-like H I. Note that the initial 
number of H I atoms, A/ki.initj is the same in both cases, 
as all the H I is high- velocity at the start of the simula- 
tion. The cool gas that accumulates along the r = axis 
after ablation, mixing, and cooling has low velocities; it 
is thus included in /?aii(£) but not /?hvc(£)- We did not 
include material that escapes from the top of the com- 
putational domain. This means that the true value of 
A/"hi(£) is larger than the value obtained from the com- 
putational domain, and so our estimates of j3(t) are upper 
limits. However, the material that escapes from the top 
of the domain has low densities and low velocities, and 
so the amounts of HVC-like material are not significantly 
affected. 

Figure [3] shows f3(i) for each of our 7 models. The black 
lines show /?aii(£)> and the colored lines show /3hvc (i). In 
all cases, f3 generally increases with time, indicating that 
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Figure [2j Continued 




50 100 150 200 250 
time (Myr) 

Figure 3. Ratio (3 of the number of H I atoms lost to the initial 
number of H I atoms, as a function of time (see Equation [TJ. The 
different line styles correspond to the different cloud models, as 
indicated in the key (thin solid line: Model A; dotted line: Model B; 
short dashed line: Model C; dot-dashed line: Model D; triple dot- 
dashed line: Model E; thick solid line: Model F; long dashed line: 
Model G). For each model, the black line shows j3 for H I with all 
velocities (/3aii), while the colored line shows ft for H I with HVC- 
like velocities (/3hvc)> respectively. The smooth thick light gray 
curves show (3(t) for spherical clouds with initial radii of 20, 150, 
and 300 pc (top to bottom) tha t lose mass at rates proportional to 
their surface areas (see Section f3. 3. 31 for details). 

H I is lost throughout the simulation. For Models A 
through E plus G, the black and colored lines are similar 
to each other (they are almost identical for Model A), in- 



dicating that most of the ablated material is "hot" (i.e., 
above 10 4 K). However, for Model F the colored line is 
clearly above the black line, indicating that some of the 
ablated material is not yet ionized, or that some of the 
ablated gas has been heated and subsequently cooled. 
We find that the latter explanation is the more impor- 
tant: the amount of radiatively cooled ablated gas that 
accumulates along the r = axis is larger in Model F 
than in the other models. 

Here we note a number of features from Figure 02 (1) 
j3 is relatively insensitive to the cloud's initial velocity, 
over a wide range of velocities (100-300 km s _1 ; compare 
Models B, C, and D). (2) /3 is also relatively insensitive to 
the cloud's initial density profile (compare Models C and 
E) and the cloud and ISM's initial densities (compare 
Models B and G). (3) A smaller cloud loses its H I con- 
tent, as a fraction of its initial mass, faster than a larger 
cloud (compare Models A, B, and F). With these various 
trends in mind, we will discuss the effects of the different 
model parameters in more detail in the following section. 

3.3. Differences between the Models - The Effects of 
Different Model Parameters 

Figure [4] compares the evolution of Models A, C, D, 
E, F, and G by plotting the temperature on a logarith- 
mic scale at various times (the evolution of Model B is 
shown in Figure [2j). Note that for Models C through 
G, the four panels correspond to similar stages in these 
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Figure 4. Comparison of the evolution of the cloud in Models A, C, D, E, F, and G (cf. Figure 0. Cross-sections through the clouds 
showing the logarithm of the temperature are plotted for Model A (top left), Model C (top right), Model D (middle left), Model F (middle 
right), Model E (bottom left), and Model G (bottom right). For Model A, the panels correspond to t = 15, 30, 45, and 60 Myr, respectively, 
from left to right. For Models C, D, E, and G, the panels correspond to t = 30, 60, 90, and 120 Myr, and for Model F they correspond to 
t = 60, 120, 180, and 240 Myr. The same color scale is used for all plots (see the color bar in the extreme top right). Models C, E, and G 
have the same size domain as Model B (Figure [2}. The Model A domain has a smaller height and width than the Model B domain, while 
the Model F domain has a larger height and width. The Model D domain has the same width as but a larger height than the Model B 
domain. 



clouds' evolution, while the first Model A panel corre- 
sponds to a similar stage of evolution as the final panels 
for the other models (see Section l3.3.3j) . The computa- 
tional domain for Model D has a larger height than those 
for Models B, C, E, and G, because of the higher initial 
velocity of the cloud. In simulations with higher initial 
cloud velocities (Models C, D, and E), the location of 
the cloud shifts upward further than in the simulations 
with lower cloud velocities (Models B, F, and G). This 
upward shift is because, in our simulations, the cloud is 
initially stationary while the ISM flows upward, pushing 
the cloud upward (see Section [2]). However, despite the 



low velocity in Model A, the cloud in this simulation is 
still shifted upward significantly, due to its low inertia 
relative to the ISM. 

3.3.1. The Effect of the Cloud Velocity (Models B, C, and 

D) 

The initial conditions for the cloud and the ISM are the 
same in Models B, C, and D, apart from the initial ve- 
locity of the cloud: —100, —150, and —300 km s _1 in the 
observer's frame, respectively. Because the sound speed 
of the ISM (T = 10 6 K) is -150 km s"\ these veloc- 
ities correspond to the subsonic, transonic, and super- 
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sonic regimes, respectively. We would therefore expect 
a bow shock to develop in Model D (supersonic case); 
Figure |4] shows that a bow shock does indeed develop at 
early times in this model, and persists until the end of 
the simulation. 

Apart from the formation of a bow shock in Model D, 
Models C and D both evolve with similar hydrodynam- 
ical processes that were seen in Model B: Bernoulli's ef- 
fect, ablation of the cloud's material due to shear insta- 
bilities, mixing of the ablated gas with the hot ISM, and 
cooling of the mixed gas. However, the shear instabilities 
grow more rapidly the larger the velocity differenc e be- 
tween the cloud and the ISM f Chandrasekharl fl 96 1 . Sec- 
tion 101). The faster-growing instabilities clearly affect 
the evolution of the cloud: a faster cloud disrupts more 
violently than a slower cloud. The temperature plots 
show that the amount of cool cloud material (shown in 
white) is smaller for a faster cloud at a given time: com- 
pare the 3rd, 5th, 7th, and 9th panels in the Figure [2] 
temperature plot (Model B) with the Model C and D 
panels in Figure |4] (these sets of panels correspond to the 
same times: t = 30, 60, 90, and 120 Myr). These differ- 
ences in the severity of the cloud disruption due to shear 
instabilities lead to different morphologies for the clouds 
with different velocities. 

Although different cloud velocities lead to different 
cloud morphologies, we find that the faster-growing in- 
stabilities do not necessarily lead to material being ab- 
lated and/or ionized at larger rates for the faster clouds. 
Figure [3] shows that the rates at which the Model B, C, 
and D clouds lose their H I to ablation and/or ioniza- 
tion, offset by material that has cooled, are rather simi- 
lar. Note that in Model D, the bow-shock that forms in 
front of the cloud helps protect the cloud from ablation, 
because the shocked ISM has a lower velocity than the 
unshocked ISM. Note also that in Model D, the mixing 
of the ablated gas with the hot ISM is also different from 
the mixing in Models B and C: the fast-moving ISM con- 
strains the ablated material so that it remains close to 
the cloud. This gas mixes and cools as it moves along the 
cloud periphery. The mixed gas is so closely constrained 
to the edge of the cloud that it cannot be clearly seen in 
the temperature plots (Figured]), but high ions are abun- 
dant along the periphery of the cloud where the mixed 
gas exists (see Figure [5]). 

3.3.2. The Effect of the Cloud Density Profile (Models C 
and E) and the Cloud Density (Models B and G) 

Models C and E have the same initial conditions, apart 
from the initial density profile of the cloud: in Model C, 
the cloud density decreases smoothly at the edge of the 
cloud until it equals the ambient density, whereas in 
Model E the cloud has a uniform density and a sharp 
edge (see Figure [T]). Models E and C allow us to com- 
pare the results for a uniform, sharp-edged density profile 
with those for a more realistic density profile. Note that 
similar comparisons have been made in previous works 
(jVieser fc Henslerl 120075 iHeitsch fc Putmanll2009D . 

The first two panels for Model C in Figured] show that 
a tail of low-temperature gas propagates upward imme- 
diately behind the cloud, whereas Model E has not devel- 
oped such a tail at early times. These tails are composed 
of the low-density gas that was initially at the edge of 
cloud in Model C, due to the smooth transition between 




r(pc) 

Figure 5. O VI ion fraction for Model D at t = 30 Myr (the same 
time as the first Model D panel in Figure [4|). 

the cloud density and the ISM density at the cloud edge. 
Similar tails are also seen in the other models at early 
times, except in Model E. The low-density gas at the 
edge of the Model C cloud is ablated via shear instabil- 
ities more quickly than the high-density gas at the edge 
of the Model E cloud. This ablated material accumu- 
lates near the r = axis at early times after mixing 
with the ISM and radiatively cooling, and moves upward 
along this axis, forming the tail. Cool gas does eventually 
start accumulating along the r = axis in Model E, but 
at later times (see the third Model E panel in Figure H] 
corresponding to t = 90 Myr). 

Models B and G have the same initial conditions, ex- 
cept for the cloud and ISM densities being an order of 
magnitude lower in Model G. The Model G cloud gen- 
erally evolves in a similar fashion to the Model B cloud; 
this is because the ratio of densities across the cloud-ISM 
interface are the same in both models, and so the growth 
of perturbations due to the shear instabil ity should oc- 
cur at the same rate (|Chandrasekhar|[l96lL Section 101). 
However, the lower density in Model G reduces the cool- 
ing rates, so the mixed gas cools more slowly. As a result, 
very little of the ablated gas that accumulates along the 
r = axis in Model G is cool. 

As mentioned at the end of Section 13. 2\ the detailed 
differences in the hydrodynamical evolution of the cloud 
between Models C and E as well as between Models B 
and G do not affect j3 significantly (Figure [3]), indicating 
that the rate at which the cloud loses its H I does not 
strongly depend on either the initial density profile or 
the initial density of the cloud, provided that the density 
contrast between the cloud and the ISM is the same. 

3.3.3. The Effect of the Cloud Size (Models A, B, and F) 

Models A, B, and F have the same initial conditions, 
apart from the size of the cloud; in these three models, 
the radii of the initial clouds are ^20, ^150 and ^300 pc, 
respectively (see Figure [1]). As these clouds all have the 
same initial velocity ( — 100 km s _1 ), they should all be 
subject to the same physical processes. 
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As noted at the end of Section 13.21 a smaller cloud 
loses its H I content more rapidly than a larger cloud. 
In order to better understand this trend, we consider the 
simple case of a uniform spherical cloud losing mass at 
a rate proportional to its surface area, i.e., dM(t)/dt oc 
— 47rr 2 (t). In this case, the radius decreases linearly with 
time, r(t) = ro — kt, where ro is the initial radius of the 
cloud and k is the rate at which the radius decreases. 
Because the cloud in this simple model is uniform and 
spherical, J\fm(t) oc r 3 (£), and so /3(t) = 1 — (1 - kt/ro) 3 
(from Equation [1} . 

The smooth thick light gray curves in Figure [3] show 
f3(t) calculated according to this simple model with k = 
0.1 pc Myr -1 and r = 20, 150, and 300 pc (top to bot- 
tom). This simple model represents the loss of H I from 
the cloud due to all processes, and so should be compared 
with the colored (HVC-like H I) curves; as previously 
noted, these curves represent H I lost to ablation and to 
ionization. Note that this model is not an accurate phys- 
ical model - Figures [2] and |4] show that the cloud does 
not remain spherical during its evolution. As a result, we 
did not fit the thick light gray curves in Figure [3] to the 
corresponding curves derived from the hydrodynamical 
simulations. Nevertheless, this simple model can provide 
some insight into the relative behavior of Models A, B, 
and F. 

The curve derived from this simple model for ro = 
20 pc is in reasonably good agreement with the Model A 
curve. For Models B and F, the simple model overes- 
timates the mass-loss at earlier times, and for Model F 
it slightly underestimates the mass- loss at later times. 
Despite these shortcomings, this simple model indicates 
that a major reason that a smaller cloud loses its mass 
more rapidly is because it has a larger surface area rela- 
tive to its mass. 

If the clouds evolve according to this simple model, 
with the same value of k for all clouds, then the Model A 
cloud will be at a similar phase in its evolution (i.e., at 
the same value of r(t)/ro) at t = 16 Myr as the Model B 
cloud at t = 120 Myr and as the Model F cloud at t = 
240 Myr. We find that the morphology of the cloud in all 
three models varies in a similar fashion according to this 
timescale; i.e., the Model A, B, and F clouds have similar 
shapes at t = 16, t = 120, and t = 240 Myr, respectively. 
To see this, compare the first Model A panel in Figure |4] 
(t = 15 Myr), the final temperature panel in Figure [2] 
(Model B at t = 120 Myr), and the final Model F panel 
in Figure!! (£ = 240 Myr). Therefore, the Model A cloud 
after t = 16 Myr (second through fourth Model A panels 
in Figure |4|) is at a later phase in its evolution than the 
Model B cloud at t = 120 Myr and the Model F cloud at 
t 240 Myr. 

At a given time t, although a smaller cloud will have 
lost more of its H I mass relative to its initial mass than 
a larger cloud, the above simple model predicts that the 
larger cloud will have lost more mass. Our simulations 
bear out this prediction - the total mass lost from the 
Model F cloud is larger than the mass lost from the other, 
smaller model clouds. Furthermore, if the ablated mate- 
rial were to mix and cool at the same rate in all models, 
then the fraction of the ablated gas that is cool would 
be the same in all models, regardless of the cloud's ini- 
tial size. This would result in the ratio of Paw to Puvc 



being the same in all models. However, Model F, which 
has the largest cloud, yields larger ratios than the other 
models; i.e., the difference between the black and colored 
lines for Model F in Figure [3] is larger than for the other 
models. This probably indicates that more radiatively 
cooling takes place in the ablated material in Model F 
than in the models with smaller clouds. 

3.4. The Fate of High- Velocity Clouds 
Hei tsch fc Putmanl (|2009f ) modeled HVCs of various 
sizes (initial H I mass, M H i,i n it = 10 3 1 -10 4 6 M©), and 
found that HVCs with M H 'i,init < 10 4 5 M© will lose 
all their H I after traveling <10 kpc. Thus, smaller 
HVCs are unlikely to reach the disk as neutral hydro- 
gen; nor will larger HVC complexes, if they are in fact 
composed of small cloudlets. However, it is possible 
that the HVC material could still reach the disk in the 
form of warm ionized material (Heitsc h fc Putm an 2009; 
IShull et al.ll2009HBland-Hawthornll2009D . 

Our suite of models include one cloud (Model F) with 
^Hi,init = 10 5,5 M© 3 w hich is larger than the ma sses of 
the clouds simulated by iHeitsch fc Putmanl (|2009[ ) . The 
Model F cloud has /3hvc ~ 0.3 at the end of the sim- 
ulation (t = 240 Myr); i.e., ^70% of the cloud's initial 
high- velocity H I mass remains at T < 10 4 K at this time, 
although it is possible that some of this high- velocity ma- 
terial has broken off the main cloud. The cloud will travel 
much further before completely dissipating - the value of 
Phvc from the end of the Model A simulation and the 
scaling discussed in Section 13.3.31 imply that ^30% of 
the Model F cloud's initial high-velocity H I will remain 
below 10 4 K at t ~ 900 Myr. 

As well as our simulating a more massive 
cloud, there are other differences between our and 
IHeitsch fc PutmanT s simulations: we used a 2D geome- 
try and a cooling curve calcu lated with solar abundances, 
whereas Hei tsch fc Putmanl carried out 3D simulations 
with 1/10-solar abundances. Both of these differences 
would tend to stabilize our clouds against disruption 
(Section [5]), leading to longer-lived clouds in our simu- 
latio ns. However, com paring similar-sized clouds in our 
and IHeitsch fc Putmanr s simulations indicates that the 
cloud lifetimes agree within a factor of ~3-5. Therefore, 
even taking this into account, our Model F simulations 
indicate that very large clouds (>10 5-5 M©) will live for 
at least a few hundred megayears, and travel a few tens 
of kiloparsecs (assuming a speed of ~100 km s _1 ). Such 
large cloud masses are not implausible: e.g., masses of 
>2 x 10 6 M©, -10 7 M©, a nd ~10 7 has been mea- 
sured for the Smith Cloud ([N ichols & Bland-Haw thorn 
120091) . Complex C (Wakke feTaLl 120071 : iThom et al. 
2008), and Complex H (Lockman 2003), respectively 
(although the complexes are not single clouds). The 
largest HVCs may therefore survive as far as the 
Galactic disk at least partially as neutral hydrogen. 

3.5. High Ions 

In our simulations, the gas that ablates from the cold 
cloud and mixes with the hot ISM is rich in high ions 
(C IV, N V, and O VI). In this section we investigate the 

4 Note that the total cloud masses in Table [T] have to be divided 
by 1.4 to give the H I masses. 
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Figure 6. Masses of C IV (top), N V (middle), and O VI (bottom) 
as a function of time from our various model simulations. The black 
lines show the total mass of each ion in the simulational domain, 
and the colored lines show the masses of the ions moving with 
HVC-like velocities. The same line styles are used for the various 
models as in Figure 03 

properties of these high ions in more detail. First, we 
estimate the quantities of high ions that are produced via 
ablation and ionization of the cloud material, and then 
we discuss the variation of column density with radius 
for each high ion; these column densities can be directly 
compared with observations. 

Figure [6] shows the masses of the high ions in our sim- 
ulational domains as functions of time. In each case we 
calculate two values: the total mass of a given ion in the 
simulational domain, regardless of velocity (black lines), 
and the mass of a given ion moving with HVC-like ve- 
locities (colored lines). When calculating the total mass 
of a given ion, integrated over all velocities, we subtract 
off the mass of that ion that was contained in the hot 
ambient ISM at t = 0. Our values do not include the 



ions that have escaped from the domain. This means 
that the total masses of the high ions are actually lower 
limits. However, as with the amount of HVC-like H I 
(Section I3.2p . the masses of the HVC-like ions are not 
significantly affected by our neglecting material that has 
flowed off the domain, because such material has a low 
density and a low velocity in the observer's frame. Note 
from Figure [6] that the masses of the high ions with HVC- 
like velocities are smaller than the corresponding total 
masses integrated over all velocities. This is because the 
gas that ablates from the cloud decelerates as it mixes 
with the ambient gas. Eventually, the mixed gas slows 
to halo-like velocities, but still contains high ions. 

In Models B through G, the masses of the high ions 
generally increase with time, although there are fluctu- 
ations on small timescales. Because the high ions are 
abundant in the material ablated from the cloud, the 
time evolution of their masses is similar to that of the lost 
H I in two distinct ways. Firstly, similar amounts of high 
ions are produced regardless of the cloud's initial velocity 
(compare Models B, C, and D; cf. Section [3. 3. ip or den - 
sity profile (compare Models C and E; cf. Section 13. 3. 2|) . 
Typically, the amounts of high ions produced in Mod- 
els B through E agree within a factor of 2. Secondly, 
more high ions are produced from the larger cloud than 
from the smaller cloud (compare Models A, B, and F; cf. 
Section [3X3]). The density of the Model G cloud is 1/10 
that of the Model B cloud, and so the masses of the high 
ions produced in Model G are commensurately lower. 
Unlike the other models, Model A shows a decrease in 
the mass of high ions at later times. This occurs because 
Model A is sampling a much later phase of the cloud 
evolution (see Section 13. 3. 3p . such that by the end of the 
Model A simulation, the cloud is mostly destroyed and 
so provides fresh cool gas for mixing with the ISM more 
slowly. 

To investigate the spatial distributions of the high ions, 
we calculated column densities for sightlines running ver- 
tically through the Model B domain. Figure [71 shows H I, 
C IV, N V, and O VI column densities as functions of ra- 
dius for t = 30, 60, 90, and 120 Myr. The footprints of 
the high ion distributions increase with time, especially 
for the ions with slow velocities. The column densities of 
the high ions generally increase with time while the H I 
column density decreases, because the cool H I-rich gas 
is ionized and converted to high ion-bearing gas. Most 
of the material containing high ions with HVC-like ve- 
locities was recently ablated from the cloud and ionized 
to this level, which is why the footprints of the HVC- 
like high ions (gray lines in Figure [7]) are similar to the 
footprint of the H I cloud. 

The column densities of H I calculated throughout the 
Model B simulation are above the current 21-cm detec- 
tion limit (a few times 10 18 cm -2 ) over most of the cloud, 
so if the Model B cloud were a real cloud, it could be iden- 
tified as an HVC throughout the whole simulation period. 
We find that simulations with similar clouds (i.e., Mod- 
els C, D, and E) have similar H I column densities to 
Model B. The H I column densities for larger (Model F), 
smaller (Model A), or less dense (Model G) clouds than 
the Model B cloud vary according to their size or den- 
sity. Even the smaller column densities (from Models A 
and G) are above the current detection limit over much 
of the clouds up to the ends of the simulations. 
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Figure 7. Column density as a function of radius for H I (top left), C IV (top right), N V (bottom left), and O VI (bottom right), from 
Model B. In each plot, the panels correspond to t = 30, 60, 90, and 120 Myr (top to bottom). The black lines are the column densities 
integrated over all velocities, and the gray lines for the column densities for material with HVC-like velocities. Note that for H I (top left), 
these two lines are almost identical. 



At the edges of the clouds, however, the H I can become 
undetectable, while the high ions remain detectable. For 
example, in Model B at t = 120 Myr, the H I col- 
umn density falls below the 21-cm detection limit be- 
yond r « 300 pc, but there are substantial amounts of 
O VI (up to ~10 13 5 cm -2 ) out to larger radii (see Fig- 
ure [7j). Also, it should be noted that the column densi- 
ties plotted in Figure [7] are for vertical sightlines through 
our model domains. The high- velocity high ions reside 
mainly behind the main body of the cloud, in material 
that has ablated from the cloud, mixed with the ambient 
medium, and fallen behind the cloud. Therefore, diag- 
onal sightlines through the model domain could inter- 
sect significant column densities of high ions, while miss- 
ing most of the H I. Our simulation results could there- 
fore partially explain the presence of highly ionized high- 
velocity gas on sightlines withou t corresponding high- 
velocity H I ([Sembach et al.ll2003[ ). Our results offer only 



a partial explanation because they account only for sight- 
lines with high- velocity high ions that have high- velocity 
H I nearby; some of the high-velocity O VI detections 
in iSembach et al.l (j2003[ ) are on sightlines without any 
high- velocity H I within several degrees. 

We noted above that similar amounts of high ions are 
produced in Models C, D, and E as in Model B (Fig- 
ure [6]) . The top four panels in Figure [8] show that the 
O VI column densities from Models C, D, and E are sim- 
ilar to those from Model B at t = 120 Myr, although the 
footprints vary from model to model. In contrast, the 
Model G O VI column densities (5th panel in Figure [8]) 
are about an order of magnitude smaller than those from 
Model B, which is commensurate with the density differ- 
ence between these two models. 

As was noted in Section I3.3.3| Model A is mostly in 
a later evolutionary phase than the other models. From 
Figure [9j we see that the O VI column densities in this 
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Figure 8. O VI column density as a function of radius from 
Models B, C, D, E, G, and F (top to bottom). All plots are for 
t = 120 Myr, except for the Model F plot, which is for t = 240 Myr. 
Note the larger r axis range for Model F. As in Figure [7] the black 
lines show the total column densities, regardless of velocity, and 
the gray lines show the column densities of HVC-like material. 

model peak around t = 30 Myr and subsequently de- 
crease (see also the bottom panel in Figure [6]) . When we 
compare Models A, B, and F at similar stages of their 
evolution (t = 15, 120, and 240 Myr, respectively; see 
Section I5.3.3|) . we find that the O VI column densities 
with HVC-like velocities from Models A and F are factors 
of ~5 smaller and ~1.5 larger than those from Model B, 
respectively (-2 x 10 13 , -1 x 10 14 , and -1.5 x 10 14 cm -2 
for Models A, B, and F, respectively). The ratio of the 
O VI column densities from Models A, B, and F (approx- 
imately 0.2:1:1.5) is similar to but not equal to the ratio 
of initial radii, 0.13:1:2. 

In the following section, we compare our column den- 
sity predictions with observations of Complex C. In par- 
ticular, we will consider the ratios of the high ions to 
H I and to each other. Such ratios predicted by different 
models are frequently compared with the observed ratios 
in order to identify the physical process (es) by which the 
high ions are produced. 

4. COMPARISON WITH OBSERVATIONS OF COMPLEX C 

Among the several HVC complexes, Complex C has 
been studied most extensively, perhaps because of 
its location in the northern h emisphere and its large 
size on the sky (—1800 deg 2 ; Wak ker fc van Woerdenl 
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Figure 9. O VI column density as a function of radius from 
Model A at t = 15, 30, 45, and 60 Myr (from top to bottom; 
these are the same times as the Model A panels in Figure 3]) As 
in Figure the black lines show the total column densities, re- 
gardless of velocity, and the gray lines show the column densities 
of HVC-like material. 

119911 ). Its distance is r eason ably well constrained (10 ± 
2.5 kpc; iThom et all 120081), and it is of r ather low 
metallicity (-0.13 solar; ICollins et all 12007ft . Obser- 
vations of high-velocity high ions along several sight 
lines passing through the Complex C region ha v e bee n 
repo rted bv ISembach et all (|2003ft . iFox et al.l ([2004( 1. 
and ICollins et al.l ([2007D . In addition, Complex C is 
thought to be coll iding with the Milky Way's thick disk 
(jTripp et al.ll2003ft . In this section, we will compare our 
model predictions with observations of Complex C in sev- 
eral different ways. First, in Section |47TJ we compare the 
ion column densities, iV, predicte d by our models with 
those measured for Complex C ( Sembac h et al.l 120031: 
IFox et all 120041 : ICollins et all 12007ft . [Sembach et al. 
(|2003ft looked for but did not find a correlation between 
the Complex C O VI and H I column densities. We dis- 
cuss this lack of correlation in terms of our model predic- 
tions in Section l4~2l Finally, in Section [4731 we calculate 
the ratios of column densities averaged over the whole 
of the model cloud, and we compare these ratios with 
observations in Sections 14.41 and 14.51 

4.1. Ion Column Densities 

ISembach et al.l (|2003l ) searched for high-velocity O VI 
in the survey of high-latitude Far Ult raviolet Spec- 
trosc opic Explorer (FUSE) observations ( Wa kker et al.l 
12003ft . They detected a total of 84 high- velocity O VI 
absorption features on 59 out of the 102 sightlines in 
the survey, spread over the high-latitude sky. The me- 
dian O VI column density of these 84 detected features is 
9.3 x 10 13 cm -2 , although there is significant sightline-to- 
sightline variation (b etween ~10 13 and ~3 x 10 14 cm -2 ; 
ISembach et al.ll2003[ ). For the 9 sightlines toward Com- 
plex C with detections, the median high-velocity O VI 
column density i s 7.6 x 10 13 cm" 2 (range: (4.7-16.6) x 
10 13 cm" 2 ; ISembach et~aTll2003ft . In Model B, the col- 
umn densities of HVC-like O VI approach the median of 
the observed column densities at some times for some 
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Figure 10. Observed and predicted O VI versus H I column den- 
sities (plotted on logarithmic scales). The filled stars sho w the 
Complex C observations (Table 8 from Sembach et al. 2003). The 
other symbols show predictions from our simulations for material 
with HVC-like velocities. The colors and symbols respectively de- 
note the different models and different times (e.g., orange diagonal 
crosses denote Model F at t = 240 Myr). For each model and 
for each time, column densities were extracted for multiple verti- 
cal sightlines; the results for each sightline are plotted. Note that 
not all the times in the key were used for each model. We consid- 
ered column densities of H I and O VI above 10 18 and 10 12 cm -2 , 
respectively. 

impact parameters, but are generally lower than the 
observed median. Furthermore, along some sight lines 
toward Complex C, there are measurements of C IV 
and N V column densities ([Fox et al. 2004; iCollins et al.l 
l2007t ). The observed high- velocity C IV column density 
ranges from <6.9 x 10 12 to 6.5 x 10 13 cm -2 while that 
ofNV ranges from <6.9 x 10 12 to <3.1 x 10 13 cm -2 
([Collins et al.ll2007[ ). For certain choices of age and im- 
pact parameter, our Model B results match the column 
densities at the higher end of the observed ranges (e.g., 
r - 150 pc at t = 120 Myr for C IV and r - pc 
at t = 120 Myr for N V) and at the lower end (e.g., 
r - 150 pc at t = 60 Myr for C IV and r > 100 pc at 
t = 60 Myr for N V). 

4.2. N(0 VI) versus N(H I) 

The filled stars in Figure [TO] show N(0 VI) 
against 7V(H I) for 9 sig htlines through Com- 
plex C (ISembach et all I2003L Table 8; we followed 
iSembach et al.l and did not include the H I components 
in brackets in their table). As noted bv ISembach et all 
there is no correlation between N(0 VI) and iV(H I); 
7V(0 VI) varies by ^0.5 dex over Complex C while 
7V(H I) varies by «1.5 dex. The other symbols in 
Figure [10] show the predictions from our various models, 
extracted for multiple vertical sightlines through our 
model domains at several different times. For any given 
model at a given time, the predictions for different 
sightlines exhibit a similar trend to the observations; 
i.e., N(0 VI) is approximately constant over a large 
range of N(H I). The trend in the predictions is because, 
in our simulations, the H I column density decreases 
significantly toward larger radii, whereas the O VI 
column densities remain more or less constant out to 
larger radii (see Figure [7j). Not only do the predictions 
exhibit the same trend as the observations, but also 
some of the predicted values are in good agreement with 
the observations. However, not all of the observations 
can be matched by a single model at a specific time. 



For example, 5 of the 9 observed data points (the four 
with the largest values of N(0 VI) and the one with 
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similar to values from Model C at t = 120 Myr (cyan 
squares) and from Model F at t = 240 Myr (orange 
diagonal crosses). The remaining data points are similar 
to predicted values at other times and/or from other 
models. 

The predicted values from Models A and G (small and 
low-density clouds, respectively) have H I and O VI col- 
umn densities that are smaller than the observed values 
for Complex C, regardless of simulation time and sight- 
line. However, it is possible that Complex C is composed 
of many small or low-density clouds, each of which indi- 
vidually resembles a Model A or G cloud. If this is the 
case, the total column densities would be obtained by 
multiplying the Model A or G predictions (which are for 
a single cloud) by the number of clouds along the line of 
sight, resulting in a shift toward larger values of 7V(H I) 
and N(0 VI). This means that the predictions for mul- 
tiple Model A or G clouds would belong upwards and to 
the right of the Model A or G points in Figure [TUJ and 
therefore could become consistent with the observations. 
Note that, in order to carry out this shift, we must as- 
sume that each cloud contributes the same amount of 
O VI and H I at the same period in time. 

The consistency between the observed and predicted 
values of N(0 VI) against iV(H I) supports the idea that 
Complex C has interacted with the hot ISM to produce 
O VI. In particular, this consistency confirms that the 
loss of cloud material due to shear instabilities and its 
subsequent mixing with the hot ISM, subject to radiative 
cooling, is a viable physical process for the production of 
O VI. However, the timescale for Complex C's interac- 
tion with the hot ISM is uncertain, as it is dependent on 
the size and density of the constituent clouds - a set of 
smaller clouds could produce similar values of 7V(0 VI) 
and 7V(H I) to a single larger cloud that has been evolv- 
ing for a longer time. 

In the following subsections, we will look at the evo- 
lution of the column density ratios for each model, and 
compare them with the Complex C observations. This 
will enable us to place a rough constraint on the timescale 
for Complex C's interaction with the hot ISM. 

4.3. Column Density Ratios 

For comparison with the observed column density ra- 
tios, we calculated ratios of column densities that had 
been averaged over the model clouds. We considered 
only sightlines for which the column density was above 
some cut-off, iV cu t, so that material not significantly af- 
fected by the cloud-ISM interaction is not included in the 
averages. For a given ion, the average column density, 
N(t), at time t is given by 



N(t) 



J N(t, r)a(t, r)r dr 
J a(t, r)r dr ' 



(2) 



where N(t,r) is the column density for a sightline at 
radius r, and 



r i if 

a = { ol 



N > N cut ; 
otherwise. 



(3) 
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For H I, N cut = 10 16 cm- 2 for Models B through F 
and 10 15 cm -2 for Models A and G. For the high ions, 
^cut = 10 11 cm -2 for Models B through F and 10 10 cm -2 
for Models A and G. We found that the choices of 7V cut 
did not significantly affect the average column densities 
and the resulting ratios We did not subtract off a back- 
ground column density due to the hot ambient medium: 
because the high ion fractions are small in this medium, 
the background column density is negligible. 

The evolution of the column density ratios calculated 
above is presented in Figures [IT] (ion to H I) and [12] (ion 
to ion). Because we will compare the models with Com- 
plex C observations, we plot only the ratios for HVC-like 
materia l. These cur ves were calculated using solar abun- 
dances (j Allenl 1 1 973h . but s ubsolar abundances have been 
repor ted for Complex C (|Fox et al.l l2004t iCollins et"aH 
2007). We can adjust our predicted ratios by shifting 
the curves in Figures [TT] and [12] according to the ra- 
tios of the Complex C abundances to solar abundances. 
The black vertical lines in Figures fTTl and fT2l indicate the 
amounts by which the curves should be shifted downward 
to go from solar abundances to the a bundances meas ured 
along the PG 1259+593 sightline (jFox et al.l l200i see 
Table 8 in KS10 for the relevant multiplication factors). 
However, it is important to note that the abundances 
measured along the PG 1259+593 sightlin e are the low- 
est a mong; the Complex C measurements (Collins et al. 
2007), and so the shifts indicated in Figures fTTl and fT2l 
would be the most extreme cases. It is also important to 
note that the Complex C abundance measurements use 
low ions, and thus may only be sampling cloud material, 
rather than the turbulently mixed material in which the 
high ions are produced. The turbulently mixed material 
includes a contribution from the ambient medium, which 
may have a different metallicity from the cloud. Note 
that nitrogen is more depleted along the PG 1259+593 
sightline than carbon and oxygen. 

Figures [TT] and [12] also show the obser ved ratios for 
Comp lex C from iFox et al.l (|2004f ) and ICollins et all 
(2007). In some cases, these two studies give slightly 
different column densities for the same sightline, due to 
differences in their measurement methods (e.g., the de- 
termination of the continuum and the velocity range used 
for the column density integration). In such cases, we 
plot the measurements from both studies. In the fol- 
lowing two subsections, we discuss the evolution of the 
ion-to-H I and ion-to-ion ratios that we have plotted, and 
compare these ratios with the Complex C observations. 

4.4. Ion-to-H I Column Density Ratios 

The predicted ratios of high ions to H I generally in- 
crease over time, although large fluctuations occur in 
some models at earlier times (Figure [TT]). The ratios are 
largest for the smallest cloud (Model A) and smallest for 
the largest cloud (Model F) at a given time, and so an 
observed ratio can correspond to a smaller cloud that has 
evolved for a short period of time, or a larger cloud tha t 
has evolved for a longer period of time (cf. Section 14. 2j) . 

5 It is also possible to calculate averaged column density ratios 
by calculating the column density ratios for each sightline, and 
then averaging these ratios; i.e., the ratio for ions X and Y is given 
by N(X)/N(Y), as opposed to N(X)/N(Y). We found that the 
two methods gave similar results. 
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Figure 11. Comparison of predicted high-ion-to-H I column den- 
sity ratios with Complex C observations (top to bottom: C IV, N V, 
O VI to H I, respectively). The predicted ratios are calculated for 
material moving with HVC-like velocities using solar abundances 
(Allen 1973); the different line styles and colors correspond to the 
different models shown in Figures 3 and 6 The black vertical lines 
denote the amount the curves should be shifted downward to go 
from solar abundances to the Co mplex C abund ances measured 
along the PG 1259+593 sightline (jFox et alJ l20Q4 > ). The observed 
ratios for different sightlines are also plotted (note that these val- 
ues are independent of time, and so their positions along the time 
axis a re arbitrary). The black data p oints are fr om Collins 'et al.l 
( 2007), and the gray data points from Fox et al. (2004). We also 
plot the errors on the measurements, although in most cases the 
error bars are smaller than the plotting symbols. Upper limits are 
denoted by the downward-pointing arrows. 

Either interpretation would be possible without knowing 
Complex C's three-dimensional structure (i.e., whether 
Complex C is composed of many small clumps, fewer 
larger clumps, or a few large clumps mingled with many 
small clumps). 

The observed ratios vary by more than an order of mag- 
nitude from sightline to sightline. The variations in the 
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Figure 12. As Figure [TT] but showing C IV to O VI (top) and 
N V to O VI (bottom). 



predicted ratios, due to differences in the age or size of a 
cloud, are as large as the observed variations. The range 
of ratios predicted with solar abundances covers the ob- 
served range. However, when we shift the predicted ra- 
tios according to the subsolar Complex C abundances, 
the models can only match the lower observed ratios. 

Although there may be a few ways to interpret the 
observed ratios in terms of our models, we provide the 
following interpretation as a possible example. The ob- 
served N(N V)/iV(H I) ratios are almost all upper lim- 
its, and do not help to constrain the models. The 
N(CW)/N(Rl) and N(0 Vl)/iV(H I) ratios observed 
toward PG 1259+593 are similar to those predicted by 
Model A at t ~ few x 10 Myr or by Models B through E at 
t ~ 100 Myr, after the abundance shifts indicated by the 
vertical bars in Figure [TT] have been applied. Nearly all 
the other sightlines have larger observed N(C IV) /N(R I) 
and N(0 VI) /N(R I) ratios than the PG 1259+593 sight- 
line. As the predicted ratios rise with time, the model 
predictions (after the abundance shift) will match these 
other sightlines at even later times. As noted above, the 
abundance shifts indicated in Figure [TT] are the most ex- 
treme cases, as the PG 1259+593 sightline has the lowest 
abundances of the Complex C sightlines. However, even 
the solar-abundance predictions from some models (e.g., 
B through E) do not match the ratios observed along 
some sightlines until t ~ 100 Myr. These results suggest 
that the models require at least ~100 Myr of evolution to 
match the observed ratios, if the radii of the composite 
clouds were initially >150 pc (Models B through E). As 
the model clouds' initial velocities were 100-300 km s _1 
(see Table [T]), this implies that Complex C has moved 



> 10-30 kpc from its initial location to its current loca- 
tion, suggesting that Complex C likely originated as ex- 
tragalactic material. Note, however, that the measured 
metallicity of Complex C is too high for a purely e xtra- 
galactic origin (jTripp et al .1120031 : iCollins et al.ll2QQ3h and 
requires a metal-enrichment process. In preliminary cal- 
culations, we have found that mixing with the Galactic 
gas noticeably enhances the metallicity of the H I HVC 
gas, but a full discussion is beyond the scope of our cur- 
rent study. 

As noted above, there is significan t sightline-to- 
sightli ne variation in the observed ratios. iSembach et al.l 
(2003) reported a weak trend for the ratio of O VI to 
H I to increase toward lower Galactic longitude for 9 
sightlines through Complex C. Complex C is oriented di- 
agonally on the Galactic coordinate grid, running from 
low latitudes and longitudes to high latitud es and lon- 
gitudes (e.g., Wakke r fc van Woerdenl Il991h . with the 
higher-latitude, higher-longit ude region of the complex 
being higher above th e disk dThom eTaIll200l . Hence, 
the trend reported by Se mbach et al.l (j2003[ ) means that 
the ratio of O VI to H I increases toward lower Galac- 
tic latitude and lower height above the disk. The ratio 
of C IV to H I also increases toward lower latitudes, al- 
though there are fewer data points for this trend (the 
N V measurements are mostly upper limits). 

This variati on of N(Q yi)/7V (H I) with longitude was 
discussed by iTripp et all (j2003[ ) , who suggested that the 
lower- longitude, lower-latitude region of Complex C is 
closer to the disk, where the ambient medium is denser. 
They suggested that this greater density leads to a more 
vigorous interaction between the ambient medium and 
the HVC, leading to greater production of O VI and 
greater ionization of H I. Our simulations are unsuitable 
for directly examining this suggestion, as we have not 
investigated the effect of an increasing ambient density 
on high ion production. An increasing ambient density 
will lead to both greater ram pressure stripping of ma- 
terial, and faster-growing shear instabilities, because of 
the lower den sity contrast between the cloud and ambi- 
ent medium (iChandrasekharl LL961I . Section 101). Both 
of these effects should lead to greater production of high 
ions. 

An alternative explanation for the observed 
N(0 VI) /AT (HI) trend is that the region of Com- 
plex C nearer the disk is at a later stage in its evolution, 
having passed through more ambient medium, and 
so has a higher high-ion-to-H I ratio (see Figure [TT]) . 
However, the range of he ights spanned by Complex C 
(~2 kpc, using data from lThom et al.ll2008[ ) divided by 
a velocity of ~100 km s _1 corresponds to a timescale 
of only ~20 Myr. From Figure [TT1 we can see that this 
timescale is insufficient to explain the range of observed 
ratios, which vary by more than an order of magnitude 
from sightline to sightline. The trend could also be 
explained by the region of Complex C toward lower 
longitudes and latitudes being initially composed of 
smaller clouds than the higher- longitude, higher-latitude 
region, which give rise to larger high-ion-to-H I ratios 
(see Figure [TT]). However, there is no obvious reason 
for the size of the clouds of which Complex C was 
composed when it began interacting with the ISM to 
vary systematically with position. 
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Another possible explanation for the observed 
7V(0 VI)/iV(Hl) trend is that the turbulently mixed 
material toward lower longitudes and latitudes has a 
higher metallicity. Such a difference in metallicity could 
arise from the mixing of metal-poor Complex C gas 
with the metal-rich ISM. If the metallicity of the halo 
decreases with height and/or Galactocentric distance^ 
then the lower- longitude, lower-latitude region of Com- 
plex C will be interacting with higher-metallicity gas 
than the higher- longitude, higher-latitude region. As a 
result, the turbulently mixed material toward lower lon- 
gitudes and latitudes would have higher abundances and 
thus higher ratios of high ions to H I. If this specula- 
tion is correct, then variation in elemental abundances 
with latitude may be apparent from other ions. O I, for 
example, does not show an anticorrelat ion between abun - 
dance and latitude (using data from lCollins et al. 2007), 
but these abundance measurements are probably prob- 
ing undisturbed cloud material rather than turbulently 
mixed material. More measurements of elemental (as op- 
posed to ionic) abundances for Complex C are needed to 
test whether or not the abundances vary with latitude. 

4.5. Ion-to-ion Column Density Ratios 

Unlike the ion-to-H I ratios, the observed ratios of C IV 
and N V to O VI do not vary greatly from sightline to 
sight line, nor do the predicted ratios vary systematically 
with cloud size or time. The observed C IV-to-0 VI ra- 
tios are in good agreement with the predicted values; the 
agreement is even better when the predicted ratios are 
shifted according to the Complex C abundances. For 
N V to O VI, a ll but the Mrk 279 measurement from 
iFox et all ([2004 filled gray square) are upper limits, and 
so do not strongly constrain the models. If the ratio of 
gas phase nitrogen to oxygen is equal to the solar ratio, 
then our results are in good agreement with the observed 
column density ratio. However, if the ratio of gas phase 
nitrogen to oxygen is ab out 1/10 solar, as it is on the 
PG 1259+593 sight line (IFox et al.ll200l . then our re- 
sults are overabundant in N V compared to O VI by 
about a factor of 10. Note that PG 1259+593 is more 
than 10 degrees higher in Galactic latitude than Mrk 279, 
so it is not implausible that the relative abundances along 
these sight lines could differ. 

Another way to look at the ratios of C IV and 
N V to O VI is to plot them in ]V(N V)/iV(OVl)- 
N(C IV)/N(0 VI) space (Figure [13]) . Although the ra- 
tios for ions moving with HVC-like velocities predicted 
by the simulations in this paper overlap the ion ratios 
predicted by our previous numerical models of TMLs 
(KS10), the ratios from some our models are lower 
than KSlO's, and better approach the ratios observed 
along three sight lines toward Complex C. However, as 
mentioned earlier, if our results are rescaled accord- 
ing to the Complex C a bundances mea sured along the 
PG 1259+593 sightline (|Fox et al.l 12001 . our predicted 
data points are shifted by —0.97 and —0.35 along the 

6 A bundances in the d isk ISM decrease with G alactocentric dis- 
tance ( Shaver et al. 1983; Rudolph et al. 2006), but in the halo 
the situation is unclear. However, because fountains of material 
from the disk into the halo likel y dominate the hea ting of the halo, 
especially the lower halo (e.g., Henlev et al. 2010, and references 
therein), it seems not unlikely that the metallicity of the halo de- 
creases with height and Galactocentric distance. 



x and y axes, respectively. After such a shift, some of 
our models' predictions would match the observed ratios 
(note that the observed N(N V)/N(0 VI) ratios are up- 
per limits, so our model predictions would be consistent 
with the observations, even after a shift of ~1 dex to the 
left). Figure [T3l confirms that ablation of cloud material 
followed by mixing with the hot ambient gas is a viable 
physical process for the production of high- velocity high 
ions. Note that some other models shown in Figure [T3] 
are in agreement with the observed ion ratios toward 
Complex C. However, these models do not distinguish 
between high- and low- velocity ions; such a distinction is 
needed for an accurate comparison with observations of 
HVCs. 



5. NEGLECTED PHYSICAL PROCESSES AND THEIR 
EFFECTS 

Our simulations used a 2D geometry because the mem- 
ory requirements for tracking the ionization evolution of 
carbon, nitrogen, and oxygen in three dimensions are 
prohibitive. This assumed geometry prevented us from 
modeling realistic magnetic field configurations, and so 
magnetohydrodynamic effects were ignored. We also ne- 
glected thermal conduction, external heating, such as 
photoelectric heating by Galactic or extragalactic UV ra- 
diation, and photoionization. 

KS10 discuss the effect of assuming a 2D geometry 
on turbulent mixing, compared with a 3D geometry, as 
well as the effect of a magnetic field. Generally, tur- 
bulent mixing is more severe in 3D simulations than in 
2D, which would lead to greater production of high ions, 
while a magnetic field is known to suppress turbulence 
(see references in KS10). It is unclear whether the com- 
bination of a magnetic field and a 3D geometry would 
result in more or fewer high ions, compared with our 
2D simulations. In addition, our previous 3D MHD sim- 
ulations showed that the dynamical effect of the mag- 
netic field depends strongly on the orientation of the field 
(|Kwak et al.l l2QQ9). As a result, the formation and size of 
the TMLs, and thus the production of high ions, would 
also be affected by the orientation of the magnetic field. 
However, it should be noted that, for a pla ne-parallel ge- 
ometr y, the 3D MHD TML simulations of lEsquivel et al.l 
(2006) predicted similar high ion column density ratios 
as the 2D NEI simulations of KS10. Finally, it should 
be noted that the magnetic field strength is likely to be 
small in the upper halo, which is the region of interest of 
this study. 

Our simulation s did no t include thermal conduction. 
IVieser fc Henslerl (|2007bl lal) studied numerically the ef- 
fect of thermal conduction on the evaporation of a cool 
cloud embedded in a hot ambient medium. They found 
that thermal conduction suppresses shear instabilities by 
smoothing out the steep temperature and density gra- 
dients at the interface b etween the hot and cool gas 
(|Vieser fc Hen sler 2007b). However, in our simulations, 
except for Model E, there are not steep gradients at the 
initial interface between the cool cloud and the hot am- 
bient medium (see Figure [T]). Also, in our simulations, 
heat diffusion due to thermal conduction is unlikely to 
be important relativ e to heat diffusion by turbulence 
(jEsquivel et al.l [20061 : KS10). For the mixed gas in our 
simulations (T - 10 5 K, n ~ 10 cm ), the Spitzer 
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Figure 13. Observed and predicted C IV-to-O VI column density ratios against N V-to-O VI column de nsity ratios, plotted on logarithmic 
scales. The filled star near the center shows the average values for low- velocity high ions in the halo (Indebetouw Sz Shull 2004); these 
low- velocity ions will be discussed in Paper II. The three solid symbols below the sta r show the obser ved ratios for three sightlines toward 
Complex C, as indicated in the key ( Collin setaT] 120071 for Mrk 876 and Mrk 279; IFox et al.l [20041 for PG 1259+59 3). The pred i ctions 
fro m various model s are shown by the sh aded regions and the lines, as follows. Gray region: The TML models of Slavin et al. (1993) 
and [Esquivel et al. ( 2006) (the Esquivel et al. predictions with and without radiative cooling lie in the lower and upper parts of this 
region, respectively, while the ISlavin et al. I predictions are concentrated to the top and the left of this region; see Figure 7 in KS10). 
Blue regions: The TML model of KS10 (upper: NEI, lower: CIE). Purple region: This paper, specifically, the values obtained from 
each model by time-aver aging the ratios for t > 10 Myr. Yellow region: NEI calculations of a radiatively cool ing shock, for a range of 
magnetic field strengths (Dopita & Sutherland 1996). Cyan region: Radiative cooling of Galactic fountain g as (Shapiro & Beniamin 1993; 
Benjamin & Shapiro 199 3) • Green region: R adiatively cooling supernova remnant shells (Slavin & Cox 1993; Shelton 1998), and steady- 
state cooling flows (Edgar Sz Chevalier 1986 ) (the predicted ratios from these m odels are similar to ea ch other). Red region: Conductive 
heating and evaporation of spherical (Bohringer Sz Hartquist 1987) and planar (Borkowski et al. 1990) clouds (the predicted ratios from 
these two models are similar to each other). Red curve: Conductively evaporating cloud with photoionization included (Gnat et al. 2010, 
specifically, the model in their Figure llf). Different points along the curve correspond to different impact parameters; the end of the line 
at «(— 1.1, —0.7) corres ponds to a sightline throug h the cloud center. Black curves: Gas cooling radiat ively from a single t emperature: 
solid line: CIE cooling (Sutherland Sz Dopita 1993); dotted line: CIE cooling with solar abundances (Gnat & Sternberg 2007, Figure 5a); 
dashed line: NEI isochoric cooling with solar abundances (Gnat Sz Sternberg 2007, Figure 5e). As the gas cools, the ions ratios will move 
along the curves, from the lower left to the upper right. 
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100 km s and the energy 
injection scale Li n j > 1 pc ( note that L }n \ cor responds to 



thermal diffusion coefficient^ is ^sp ~ 10 24 cm 2 s 1 . 
In contras t, the turbulent diffusion coefficient ftturb ~ 
l^turb^inj (|Cho et al.ll2003l ) is >10 25 cm 2 s" 1 , where the 

turbulent velocity ^turb 
> 

the size of the largest eddies; ICho e t al. 2003). These dif- 
fusion coefficients mean that turbulent diffusion of heat 
would dominate over diffusion by thermal conduction, 
even if thermal conduction were included in our simula- 
tions. 

We included radiative cooling in our simulations (al- 
beit under the assumption of CIE), but not photoelectric 
heating due to UV irradiation of dust grains. Includ- 
ing such heating would, in principle, reduce the effect 
of radiative cooling. However, in practice, the photo- 
electric heating rate is unlikely to be important for our 
simulations of HVCs in the upper halo, because of the de- 
crease in the UV radiation field and the dust density with 
height. For example, the photoelect ric heating rate for 
the upper halo (\z\ > 3 kpc) used bv lJoung fc Mac Low! 
(|20Q6h in their ISM models would be <10 -32 erg cm" 3 
s _1 in our mixed gas, whereas the radiative cooling rate 
is ^10~ 25 erg cm" 3 s -1 (using the 1993 version of the 
iRaymond fc Smith! fl977l code). 

A more important cooling-related effect is that the 
lower metal abundances seen on some sightlines t hrough 
Complex C (jFox et al.l 120041 : iCollins et all I2007D would 
reduce the cooling rate, as the cooling is dominated by 
metal line emission. Radiative cooling tends to stabilize 
a cool cloud against disruption as it moves through hot 
ambient gas ( Vie tri et al.lll997l h and so reduced cooling 
rates would tend to accelerate the disruption of the cloud 
and so may lead to the production of more high ions than 
are currently predicted by our simulations. We defer in- 
vestigation of this effect to a later study. 

Finally, we did not consider radiative transfer, and the 
possibility that some of the ion populations may be af- 
fected by photoionization. Including photoionization is 
beyond the scope of this study. However, we note that, 
if it is important, photoionization by escaping Galactic 
radiation would tend to increase the amount of C IV rel- 
ative to O VI. 

6. SUMMARY 

We have presented the results of hydrodynamical sim- 
ulations of HVCs that self-consistently trace the non- 
equilibrium ionization evolution of carbon, nitrogen, and 
oxygen. This is the first time such simulations have 
been carried out with an initially spherical cloud mov- 
ing through hot, tenuous ambient gas. We have concen- 
trated particularly on the production of high ions (C IV, 
N V, and O VI) in the turbulent mixing layers that form 
between the cold clouds and the hot ambient medium. 
In contrast to the previous plane-parallel TML models 
(e.g., KS10), the new simulations enable us to investi- 
gate how the cold gas ablates from an initially spherical 
cloud and subsequently mixes with the hot ambient gas. 

Material is constantly being ablated from the clouds 
(Section 13. ip . This initially cold material mixes with 
the hot ambient medium, and becomes ionized, so the 



7 Derived from the Spitzer] (|19560 thermal conductivity, 6.1 x 
10- 7 (T/K) 5 / 2 erg cm" 2 s" 1 K _1 . 



clouds are constantly losing their H I content (Sec- 
tion [32]). However, despite this loss of material, we find 
that HVCs can survive for long times: at least a few hun- 
dred megayears for a cloud with initial mass ~4x 10 5 M . 
Note that clouds that are currentl y ^10 kpc from the disk 
(e.g., Complex C with z ~ 8 kpc: iThom et al.l l2008) and 
traveling toward the disk at ~100 km s _1 would require 
only ~100 Myr to reach the disk. It is therefore possible 
that the more massive regions of HVC complexes, if not 
composed of many small cloudlets, may reach the disk at 
least partially as neutral hydrogen (Section 13. 4|) . 

High ions are produced in the mixed gas, both by ion- 
ization of the initially cool cloud material and recombi- 
nations in the initially hot ambient material (Sections 13. II 
and !3.5p . The high ions initially appear while the mixed 
gas is still traveling at HVC-like velocities. The material 
stripped from the clouds slows while remaining abundant 
in high ions, and ultimately ends up with halo-like ve- 
locities (Paper II), ceasing to be recognizable as HVC 
material. 

The simulations may also explain why some high- 
velocity highly ionized gas is seen on sight lines that 
lack high- velocity H I ( Sem bach et al.l (2003). at least 
for sightlines that have high- velocity H I nearby (Sec- 
tion [33]). At the edges of the clouds, the high- velocity 
H I column density can fall below the 21-cm detection 
limit (iV(H I) < 10 18 cm -2 ), while the column densi- 
ties of high- velocity high ions remain detectable (e.g., 
A^O VI) > 10 13 cm" 2 ). In addition, the high-velocity 
high ions exist mainly in material that has ablated from 
the main body of the cloud and mixed with the ambient 
gas. This material continues to travel at HVC-like veloc- 
ities, although it has a tendency to fall behind the cloud. 
As a result, sightlines that cross the ablated material but 
not the cloud (such sight lines would cross our domain 
diagonally) would intersect large numbers of high ions 
but small numbers of neutral hydrogen atoms. 

We investigated a suite of models with a range of model 
parameters. We found that the cloud's initial velocity 
does not affect the rate of HI loss (Section 13.3. ip or 
high ion production (Section l3.5j) . The rate of H I loss 
also does not strongly depend on the initial density of 
the cloud or its density profile (Section 13. 3. 2p . However, 
the cloud's initial size does affect the rate of H I loss - 
a smaller cloud loses its H I mass relative to its initial 
mass more rapidly than a larger cloud. This difference 
is at least partially due to the fact that a smalle r clou d 
has a larger surface-area-to- volume ratio (Section [3. 3. 3p . 

Significant column densities of high ions are produced 
in our simulations. In some models, at some times, these 
column densities are large enough to be observed, and are 
in reasonable agreement with the ion column densities 
observed for Complex C (Section l4.ip . As the ion col- 
umn densities generally rise throughout our simulations, 
it is likely that the predicted column densities would be 
in better agreement with the observed values if our sim- 
ulations were continued beyond their current end times. 
Our models also compare well with the Complex C ob- 
servations in terms of the trend between N(0 VI) and 
N(R I) (Section [O]), the ion-to-H I ratios (Section ED, 
and the ion-to-ion ratios (Section [43]). We therefore con- 
clude that the ablation of cloud material and its subse- 
quent turbulent mixing with the hot ambient medium is 
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a viable mechanism for the production of high-velocity 
high ions. 

The column density ratios of high ions to H I depend 
on the initial cloud size. A smaller cloud reaches a given 
ion-to-H I ratio at an earlier time than a larger cloud. 
If Complex C were initially composed of clouds with 
r > 150 pc (which cannot be ruled out, given the lack 
of knowledge of Complex C's three-dimensional struc- 
ture), then the observed ion-to-H I ratios would imply 
an age for Complex C of >100 Myr. This in turn implies 
that Complex C has traveled > 10-30 kpc from its initial 
location, suggesting that Complex C formed from extra- 
galactic material. Furthermore, if parts of Complex C 
are composed of clouds as large as our Model F cloud 
(M ~ 4 x 10 5 Mq), these clouds may reach the disk par- 
tially in the form of neutral hydrogen, as noted above. 
HVCs could therefore provide new material to the disk 
for star formation. We also noted that there is a trend 
in Complex C for the O VI-to-H I ratio to increase as 
one gets closer to the disk, and speculated on possible 
reasons for this trend in Section 14.41 

In this paper, we have concentrated particularly on 
high- velocity ions in the Galactic halo. However, low- 
velocity high ions are also observed in the halo, and, as 
noted above, low- velocity ions are produced abundantly 
in our simulations. These low- velocity ions will be dis- 
cussed in Paper II. 

We thank Orly Gnat for providing us with the data 
from iGnat et al.) (|2010f ) plotted in our Figure [Ml We 
also thank the anonymous referee, whose comments 
have helped improve this paper. The software used in 
this work was in part developed by the DOE-supported 
AS C /Alliance Center for Astrophysical Thermonuclear 
Flashes at the University of Chicago. The simula- 
tions were performed at the Research Computing Center 
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